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Abstract 

In 60 pb^^ of data taken on the ip{3770) resonance with the CLEO-c detector, we find 8 
jx^v decay candidates that are mostly signal, containing only 1 estimated background. Using this 
statistically compelling sample, we measure a value of B{D'^ — /x+i/) = (3.5 it 1.4 it 0.6) x 10~^, 
and determine /^+ = (202 ± 41 ± 17) MeV. 
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I. INTRODUCTION 



Measuring purely leptonic decays of heavy mesons allows the determination of meson 
decay constants, which connect measured quantities, such as the BB mixing ratio, to CKM 
matrix elements. Currently, it is not possible to determine experimentally from leptonic 
B decays, so theoretical calculations of fs must be used. The most promising of these 
calculations involves lattice QCD [1-3], though there are other methods [4-8]. 

Measurements of pseudoscalar decay constants such as provide checks on these cal- 
culations and help discriminate among different models. 
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FIG. 1: The decay diagram for Z)+ jjl^v. 
The decay diagram for — > /i+i^ is shown in Fig. 1. The decay rate is given by [10] 

V{D^ ^t.) = ^flMMn.{l-^^\v,a? , (1) 

where Md+ is the mass, is the mass of the final state lepton, Vcd is a CKM matrix 
element equal to 0.224 [11], and Gp is the Fermi couphng constant. Various theoretical 
predictions of /d+ range from 190 MeV to 350 McV [1-8]. Because of helicity suppression, the 
electron mode e~^i/ has a very small rate in the Standard Model [9]. The relative widths 

are 2.64 : 1 : 2.3 x 10^'' for the r+z^, ii'^p and e+z^ final states, respectively. Unfortunately 
the mode with the largest branching fraction, r'^i', has at least two neutrinos in the final 
state and is difficult to detect. 



II. THE CLEO-C DETECTOR 

The CLEO-c detector is equipped to measure the momenta and direction of charged 
particles, identify charged hadrons, detect photons, and determine with good precision their 
directions and energies. Muons above 1.1 GeV can also be identified. The detector is almost 
cylindrically symmetric with everything but the muon detector inside a superconducting 
magnet coil run at a current that produces an almost uniform 1.0 T field. The detector 
consists of a six-layer wire drift chamber at small radius that is low mass, suitable for these 
relatively low energies. It is followed by a 47-layer drift chamber; both chambers use a gas 
mixture of 60% Hehum and 40% Propane. These two devices measure charged track three- 
momenta with excellent accuracy. The drift chamber also measures energy loss, dE/dx, that 
is used to identify charged tracks below about 0.7 GeV [12]. After the drift chamber there is 
a Ring Imaging Cherenkov Detector (RICH) [13], that identifies charged particles over most 
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of their momentum range. The RICH is surrounded by a ThaUium doped Csl crystal array 
consisting of about 8000 tapered crystals, 30 cm long and about 5x5 cm^ at the front [14] . 

III. DATA SAMPLE AND SIGNAL SELECTION 

In this study we use 60 pb~^ of CLEO-c data produced in e"'"e~ collisions and recorded 
at the tp" resonance (3.770 GeV). At this energy, the events consist of a mixture of pure 
D~^D~, D°D and three-flavor continuum events, resulting from the production of uu, dd or 
ss quark pairs. There also may be small amounts of r"'"r~ pairs and two-photon events. 

We examine all the recorded events and retain those containing at least one charged D 
candidate in the modes listed in Table I. The selection criteria are described in detail in 
what follows. We then use this sample to look for cases where we have only a single muon 
candidate whose four-momentum is consistent with a two-body D decay into a muon and a 
neutrino and no other charged tracks or excess neutral energy are present. 

All acceptable track candidates must have a helical trajectory that approaches the event 
origin within a distance of 5 mm in the azimuthal projection and 5 cm in the polar view, 
where the azimuthal projection is in the bend view of the solenoidal magnet. Each track 
must possess at least 50% of the hits expected to be on a track, and it must be within the 
fiducial volume of the drift chambers, | cos^| < 0.93, where 9 is the polar angle with respect 
to the beam direction. 

We use both charged particle ionization loss in the drift chamber (dE/dx) and RICH 
information to identify kaons and pions used to fully reconstruct D mesons. The RICH is 
used for momenta larger than 0.55 GeV. Information on the angle of detected Cherenkov 
photons is translated into a likelihood of a given photon being due to a particular particle. 
Contributions from all photons associated with a particular track are then summed to form 
an overall hkelihood denoted as Ci for each particle hypothesis. To differentiate between 
pion and kaon candidates, we use the difference: — 21og(£^) + 21og(£/^). Usually this cut 
is set at zero except for muon candidates where the difference — 21og(£^j) + 21og(£ft:) is 
required to be less than 10, to ensure a high, well understood efliciency. To utilize the 
dE/dx information we calculate as the difference between the expected ionization loss for 
a pion and the measured loss divided by the measurement error. Similarly, ax is defined in 
the same manner using the expected ionization for a kaon . 

We use both the RICH and dE/dx information for D~ meson tag candidate tracks in the 
following manner: (a) If neither the RICH nor dE/dx information is available, then the track 
is accepted as both a pion and a kaon candidate, (b) If dE/dx is available and RICH is not 
then we insist that pion candidates have PIDdE = o'l ~ '^k < 0) kaon candidates have 
PIDdE > 0. (c) If RICH information is available and dE/dx is not available, then we require 
that PIDrich = -21og(£^) + 2log{CK) < for pions and PIDrjch > for kaons. (d) If 
both dE/dx and RICH information are available, we require that {PID^e + PIDrjch) < 
for pions and {PIDdE + PIDrjch) > for kaons. 

We reconstruct tt^'s by first selecting photon candidates from energy deposits in the 
crystals not matched to charged tracks that have deposition patterns consistent with that 
expected for electromagnetic showers. Pairs of photon candidates are kinematically fit to the 
known tt" mass. We require the pull, the difference between the raw and fit mass normalized 
by its uncertainty, to be less than three for acceptable tt° candidates. 

Ks candidates are formed from a pair of charged pions which are constrained to come 
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from a single vertex. We also require that the invariant mass of the two pions be within 4.5 
times the width of the Kg mass peak, which has an r.m.s. width of 4 MeV. 

IV. RECONSTRUCTION OF CHARGED D TAGGING MODES 

Tagging modes are fully reconstructed by first evaluating the difference in the energy, 
/S.E, of the decay products with the beam energy. We then require the absolute value of this 
difference to be within 20 MeV of zero, approximately twice the r.m.s. width, and then look 
at the reconstructed D~ beam-constrained mass defined as 



mD = ^El^m-{i:%Y. (2) 

where i runs over all the final state particles. The beam-constrained mass has better resolu- 
tion then merely calculating the invariant mass of the decay products since the beam has a 
small energy spread. Besides using D~ tags and searching for ^ + u, we also use the 

charge- conjugate tags and search for A^"'^^! in the rest of this paper we will not 

mention the charge-conjugate modes explicitly, but they are always used. 

The rriD distributions for all D~ tagging modes considered in this data sample are shown 
in Fig. 2 and listed in Table I along with the numbers of signal events and background 
events within ±3 r.m.s. widths of the peak. The event numbers are determined from fits 
of the mu distributions to Gaussian signal functions plus a background shape. We fit with 
two different background parametrizations: (a) a 3'''^ order polynomial, (b) a shape function 
analogous to one first used by the ARGUS collaboration [20] which has approximately the 
correct threshold behavior at large niD] To use this function, we first fit it to the data selected 
by using AE sidebands, mode by mode, defined as 40 MeV < |A£'| < 60 MeV to fix the 
shape parameters in each mode allowing the normalization to float. For the K+7r~7r~7r°, 
Ksn~7r~7r'^ and Ks7r~7r° modes we use a single Gaussian to describe the signal whose mass 
and width are allowed to float. For the X+tt^tt" and Kstt~ modes, where we see a small 
tail on the higher mass side, we use the sum of two Gaussian's for a signal function [15]; in 
this case both the means and widths of both Gaussians are allowed to float. 



Mode 




Signal 


Background 


K+TT- 


'TT 


15173 ± 140 


583 


K+TT- 


'TT'TT° 


4082 ± 81 


1826 






2124 ±52 


251 






3975 ± 81 


1880 






3297 ± 87 


4226 


Sum 




28651 ± 207 


8765 



TABLE I: Tagging modes and numbers of signal and background events determined from the fits 
shown in Fig. 2. 

The difference between using the polynomial and ARGUS shapes in the signal yields is 
±2.2%, which we use as an estimate of the systematic error. 
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FIG. 2: Beam-constrained mass distributions for different fully reconstructed D~ decay can- 
didates in the modes: (a) D~ K~^tt~tt~, (b) D~ it~ it~ tt^ ^ (c) D~ Ks7r~ , (d) 
D~ KsTT~'K~TT~^ and (e) D~ Kstt~tt^. The solid curves show the sum of Gaussian sig- 
nal functions and S*"*^ order polynomial background functions. A single signal Gaussian is used for 
all modes except for modes (a) and (c) where the sum of two Gaussians are used. The dashed 
curves indicate the background fits. 
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Selecting those candidates within 3 r.m.s. widths of the D" mass reduces the signal 
number by 77 events giving a total of 28,574±207±629 single tag events that we use for 
further analysis. In the case of two Gaussians the wider width was used. 

V. D+ ix+v^ SELECTION CRITERIA 

Using our sample of D~ event candidates we search for events with a single additional 
charged track presumed to be a yu"^. Then we infer the existence of the neutrino by requiring 
a measured value near zero (the neutrino mass) of the missing mass squared (MM^) defined 
as 2 

MM^ = {E,,am-E^+f-(--[f^--^^+) , (3) 

where is the three-momentum of the fully reconstructed D~. 

We need to restrict the sample to candidate /J^'^i'i^ events resulting from the other D. 
Thus we wish to exclude events with more than one additional track with opposite charged 
to the tagged D, which we take to be the muon candidate, or with extra neutral energy. 
It is possible, in fact even likely, that the decay products of the tagging D~ interact in the 
detector material, mostly the EM calorimeter and spray tracks and neutral energy back into 
the rest of the detector. To evaluate the size of these contributions we use a very pure sample 
of events obtained by finding fully reconstructed D°D° events. The numbers of these events 
in various decay modes are listed in Table II, a total of 782 events. 



Mode 1 


Mode 2 


# of events 






89 






392 


if'*"7r~"7r+7r~ 


i^~7r'^7r'~7r"^ 


301 



TABLE IL Fully reconstructed D°D° events 

The number of interactions of particles with material and their consequences depend on 
the number of particles, the kind of particles and their momenta. Thus, the sum over these 
neutral D decay modes isn't quite the same as the sum over the tagging D~ decay; however, 
after accounting for the differences between the pion-nucleon and kaon-nucleon cross sections 
and the different momentum distributions of the tracks, we find that the average over these 
modes is quite similar to the D~ tagging modes for this level of statistics. 

Extra tracks do appear in these D°D° events. None of these tracks, however, approach 
the main event vertex. Requiring that good tracks are within 5 cm along the beam and 5 
mm perpendicular to the beam does not include any additional tracks from interactions in 
the material. We also reject tags with additional Ks tt+tt" candidates. 

In the D°D events, energy in the calorimeter not matched to any of the charged tracks is 
shown in Fig. 3. Figure 3(a) shows the energy of the largest shower and 3(b) shows the total. 
We accept only as extra showers those that do not match a charged track within a connected 
region. A connected region is a group of adjacent crystals with energy depositions which are 
nearest neighbors. This suppresses hadronic shower fragments which would otherwise show 
up as unmatched showers. Hadronic interactions and very energetic 7r°'s tend to produce 
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one connected region with many clusters. For further analysis we require that the largest 
unmatched shower not to be larger than 250 MeV. This requirement is (93.5±0.9)% efficient 
for signal events, estimated from the distribution of extra energies in the D°D tag sample. 
We assign an additional 4% systematic error, due to the difference in our double tag and 
single tag samples. 
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FIG. 3: Largest (a) and total extra shower (b) energies in the D°D° sample. In both cases the first 
bin is truncated; each plot has 782 total entries. 



The muon candidate is required to be within the barrel region of the detector | cos 6*1 < 
0.81; this requirement insures that the MM^ resolution is good as tracks at larger angles 
cross fewer tracking layers and consequently are measured with poorer precision. In addi- 
tion, this requirement helps reject background from the decay —>■ n~^7r°] this mode also 
gives a MM^ near zero. Requiring the muon candidate in the barrel region (the tt"'" in this 
case) avoids having the photons from this decay being lost in the transition region of the 
calorimeter between the barrel and the endcap, because the ir" direction is almost directly 
opposite the tt"*". Furthermore, the muon candidate is required not to be consistent with the 
kaon hypothesis using RICH information. Finally, we also require that the muon candidate 
deposits less than 300 MeV of energy in the calorimeter, characteristic of a minimum ioniz- 
ing particle. This requirement is very efficient for real muons, and rejects about 40% of the 
pions as determined using a sample of reconstructed D° K^n^ decays. Figure 4 shows 
the muon deposited energy in the EM calorimeter both from data on e~^e~ — > /i"^/^" and 
from GEANT simulation of the same process. The Monte Carlo and data are in excellent 
agreement for muon shower energies. We therefore use a GEANT simulation of — »• /x+z/ 
with lower energy muons to determine that the efficiency of the calorimeter energy cut is 
(98.7 ±0.2)%. 

When evaluating MM^ using Eq. (3) there are two important considerations that are not 



8 



onn 1630804-079 

oUUi — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — 1— 



> 

^ 400 



>1000 



CD 

i 500 



(a) 



J I I I I I L_UJ I I I I L 



J I I I I I I 



J I I I I I L 



T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 
(b) 




J I I I I L 



0.1 0.2 0.3 0.4 
CC Deposited Energy (GeV) 



0.5 



FIG. 4: Deposited energy in the crystal calorimeter of muons created in the process e'^e 
from (a) data and (b) Monte Carlo. 



obvious. First of all, we explicitly need to take into account the crossing angle between 
the e"^ and e~ beams. This angle is about 4 mrad, varying slightly run to run; we use this 
information and Lorentz transform all laboratory quantities to the center-of-mass. Secondly, 
we change the reconstructed D~ momenta so that they give exactly the known D~ mass; 
this changes and improves somewhat our knowledge of the D~ direction. 

The MM^ from Monte Carlo simulation is shown for our different tagging samples in 
Fig. 5. The signal is fit to a sum of two Gaussians with the wider Gaussian having about 
30% of the area independent of tagging mode. The resolution (a) is defined as 

a = /iai + (l-/i)c72, (4) 

where cJi and a2 are the individual widths of the two Gaussians and /i is the fractional area 
of the first Gaussian. The resolution is approximately 0.025 GeV^ consistent among all the 
tagging decay modes. 

Wc check our simulations using the —> Ksti^ decay. Here we choose events with 
the same requirements as used to search for /i+z/ but require one additional found Kg- The 
MM^ distribution for this final state is shown in Fig. 6 and peaks as expected at the Ks 
mass-squared of 0.25 GeV^. The resolution is measured to be 0.024±0.002 GeV^ from a 
single Gaussian fit, consistent with but slightly larger than the Monte Carlo estimate of 
0.021±0.001 GcV^. To account for the difference in resolution between data and simulations 
we scale the resolution by 14% to 0.028 GeV^ when looking for the — > ii^v^ signal. 

The MM^ distributions for our tagged events requiring no extra charged tracks besides 
the muon candidate and showers above 250 MeV as described above is shown in Fig. 7. We 
see a small signal near zero containing 8 events within a 2a interval, —0.056 GeV^ to +0.056 
GeV^. This signal is most likely due to the —>■ n'^Vfj, mode we are seeking. The large 
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FIG. 5: Monte Carlo simulation of — > /x+f^ events for different tags. The plots have been fitted 
to two Gaussians centered at zero where the second Gaussian constitutes around 30% of area. 



peak centered near 0.25 GeV^ is from the decay K"n^ that is far from our signal 

region and is expected since many Kl would escape our detector. 

Table III lists the properties of each muon candidate from the 8 events in the signal region. 
A typical event is shown in Fig. 8. 
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Tncr 


MM2 


CC energy 








(GeV2) of /i+(GeV) 






J\ TTTVTT 


0.032 


0.186 




1 Rfi n -1- 

-iOO.U 


Kg IX 


-0.019 


0.201 


n n 
u.u 


1 /in n 

-i4U.U 


KlTTT 


-0.051 


0.190 


31.9 


-252.9 + 


Kmr 


-0.004 


0.221 


0.0 


-115.2 + 


Kstttt" 


0.032 


0.164 


-0.3 


-130.6 - 




0.001 


0.245 


-11.7 


-138.9 + 




0.002 


0.204 


-8.6 


-88.6 - 




0.014 


0.208 


-8.3 


-113.0 + 



TABLE III: Muon Candidate Properties. (CC indicates the crystal calorimeter.) 



1630804-078 

n — I — \ — \ — 




FIG. 6: MM^ distribution for the decay Ksix^ from data and signal Monte-Carlo simulation 

VI. BACKGROUND EVALUATION 
A. Introduction 

There are several background sources we need to evaluate. These include background from 
other modes, background from misidentified D°D events and continuum background. 
The requirement of the muon depositing <300 MeV in the calorimeter, while about 99% 
efficient on muons, rejects only about 40% of pions as determined from the D°D° event 
sample where the pion from the X^tt^ mode was examined. In Fig. 9 we show the deposited 
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FIG. 7: MM^ using D~ tags and one additional opposite sign charged track and no extra energetic 
showers (see text). The insert shows the signal region for — > ii~^u enlarged; the ±2(7 range is 
shown between the two arrows. 



energy in the calorimeter for both kaons and pions obtained from the Ktt tag sample. 



B. Backgrounds 

There arc a few decay modes that could mimic the signal. These are listed in Table IV 
along with the background estimate we obtained by Monte Carlo generation and reconstruc- 
tion of each specific mode. The branching ratios are from the Particle Data Group except 
for the tt+tt" mode where a separate CLEO analysis gives a somewhat lower value [16]. This 
mode is the most difficult to reject because the MM^ peaks very close to zero, at 0.018 GeV^, 
well within our resolution of 0.028 GeV^. While we have insisted that the muon candidate be 
well within our acceptance, it is possible for the photons from the 7r° decay to inadvertently 
be matched to the tracks from the tagging D~ or be missed. The maximum photon energy 
of the n° from a GEANT simulation of D'^ — > tt+tt" is shown in Fig. 10. We note that at 
least one photon from the tt+tt" mode exceeds our 250 MeV calorimeter energy requirement 
and should in most cases cause such a decay to be vetoed. 

Even though the K ir^ mode gives a large peak in the MM^ spectrum near 0.25 GeV^, 
our simulation shows that only a very small amount can enter our signal region, only 0.06 
events. We have simulated backgrounds from — > t'^i/. Out of 10,000 simulated events 
with tags, we found background only when t+ — * 7r~^i/. Because of the small D^-t^ mass 
difference, the r"*" is almost at rest in the laboratory frame and thus the 7r+ has relatively 
large momentum causing the MM^ distribution to populate only the low MM^ region, even 
in this case with two missing neutrinos. The MM^ distribution is shown in Fig. 11. 

The semileptonic mode yr^yU+z/^ is similar to tt+tt" except that the tt" often carries off 
enough momentum to result in large MM^. We found no candidate background events in a 
Monte Carlo sample consisting of 50,000 tags plus a —>■ -n^jji^u decay. 
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Run: 202742 i630804-076 




Kgjt~jt+jr+ Tag 



FIG. 8: A typical fJ- T^^i event. The tag is this case is Kstt^tt^tt^ ■ The muon and 

the two oppositely charged pions forming the Ks are indicated. The tt~ is the "curler" track with 
momentum around 50 MeV. 



Mode 


B (%) 


# of events 


7r+7r° 


0.13±0.02 


0.31±0.04 




2.77±0.18 


0.06±0.05 




2.64X B{D+ 


/u+i/) 0.30±0.07 


7r°/i+i/ 


0.25 ±0.15 


negligible 


sum 




0.67±0.09 



TABLE IV: Backgrounds from specific decay modes 
C. D°D° and Continuum Backgrounds 

These backgrounds are evaluated by analyzing Monte Carlo samples corresponding to 5.2 
times the total amount of data in our possession. To normalize our Monte Carlo events to 
our data sample we used cr^o7j° = 3.5 nb and acontinuum = 14.5 nb [17]. In each sample we 
found one background event within two standard deviations of zero. These correspond to 
0.16±0.16 D°D events and 0.17±0.17 continuum events forming background. As a check on 
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FIG. 9: Deposited energy in EM calorimeter for (a) kaons, (b) pions from — > tt"* 



the continuum background we analyzed 23 pb ^ of continuum data taken a center-of-mass 
energy of 3670 MeV. We didn't find any — > //+z/ candidate events. 



D. Background Summary 

Our total background is 1.00±0.25 events. The probability of one background event fluc- 
tuating to 8 or more signal events is only 10~^, and even including the 0.25 event uncertainty 
in the background the signal has greater than five standard deviation significance. Because 
of the uncertainties in the Monte Carlo simulation we assign a 100% error to our background 
estimate: l.Oil.O events, for the purpose of evaluating the branching ratio. 



VII. BRANCHING RATIO AND DECAY CONSTANT 

Subtracting the 1.0 event background from our 8 events in the signal region, we determine 
a branching fraction using a detection efficiency for the single muon of 69.9%. This efficiency 
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FIG. 10: Maximum photon energy of the Tr° in the tt~^tt° decay from a GEANT simulation. 

includes the selection on MM^ within ±2(T limits, the tracking, the particle identification, 
the probabihty of the crystal energy being less than 300 MeV, and the probability of not 
having another unmatched shower in the event with energy greater than 250 MeV. We assign 
a relative 5.3% error on this efficiency, the components of which are shown in Table V. We 
use a 3% systematic error on track finding found using the double tagged events and we 
estimate the error on the particle identification cut to be 1% from studies of decays in 
higher beam energy data. The error on the minimum ionization cut on the muon candidate 
in the calorimeter is 0.2% and discussed in detail in section V. A 4% dominantly systematic 
error due to rejection of events with excess shower energy is assigned to the efficiency of this 
cut determined by using the D°D sample and also discussed in section V. 

Systematic error (%) 

MC statistics 0.8 
Track finding 3 
PID cut 1 
Minimum ionization cut 1 
Extra showers cut 4 

Total 5^3 



TABLE V: Systematic errors on the D'^ — > efficiency. 
To compute the branching ratio we use 7.0±2.8 signal events divided by 69.9% and the 
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FIG. 11: Missing Mass squared distribution for — t'^u and r+ ir^u. 

28574 tags. No other efficiencies enter. The systematic error on the branching fraction 
arises from the 5.3% systematic error on the efficiency, a 2.2% systematic error in the number 
of D~ tags and a 15.4% systematic error on the background. The total systematic error, 
evaluated by adding these contribution in quadrature, is 16.4%. Our result for the branching 
fraction is 

B{D+ i^+u^) = (3.5 ± 1.4 ± 0.6) X 10"^ . (5) 

The decay constant f£)+ is then obtained from Eq. (1) using 1.04 ps as the D'^ lifetime 
and 0.224 as \Vcd\ [H]- Our final result is 



/^+ = (202 ± 41 ± 17) MeV 



(6) 



VIII. CONCLUSIONS 

There have been several experimental studies of D meson decay constants. The Mark 111 
group published an upper limit of B{D~^ — > n^i'^) < 7.2 x 10""^, which leads to an upper limit 
on the decay constant fD+ < 290 MeV at 90% confidence level based on 9.3 pb~^ of data 
taken on the ijj" [18]. BES claimed the observation of one event at a center-of-mass energy 
of 4.03 GcV with a branching ratio of (0.08lo;o5)% [19]. Recently, using 33 pb~^ of ^" data 
they presented 3 event candidates and with an estimated background of 0.33 events where 
neither tt+tt", or r'^u were mentioned as a possible background modes, nor was continuum 
background considered [21]. Here they find a branching ratio of (0.122jlo;o53='= 0-010)%) ^ 
corresponding value of //)+ = (37llJi9±25) MeV. Our value is considerably smaller, though 
compatible with their large error. 
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Our analysis shows the first statistically significant signal for £)+ — > /I'^iy. The branching 
fraction is 

B{D+ //+z/) = (3.5 ± 1.4 ± 0.6) x 10"^ , (7) 
and the decay constant is 

/d+ = (202±41±17) MeV . (8) 

Our result for //)+ , at the current level of precision, is consistent with predictions of lattice 
QCD and models listed in Table VI. 

Model fD+ (MeV) fi,+/fD+ 



Lattice QCD (Fcrmilab and MILC) [2] 225l}J ± 21 1.17 ± 0.06 ± 0.06 

1.04 
1.02 



Quenched Lattice QCD (UKQCD) [3] 210 ± 10^^^ 1.13 ± 0.02lg:g^ 



QCD Spectral Sum Rules [5] 203 ± 20 1.15 ± 0.04 

QCD Sum Rules [6] 195 ± 20 

Relativistic Quark Model [7] 243 ±25 1.10 
Potential Model [4] 238 1.01 

Isospin Mass Splittings [8] 262 ± 29 

TABLE VL Theoretical predictions of 7^3+ and f^+/fjj+ 

The models generally predict /^+ to be 10-15% larger than fu+. CLEO previously mea- 
sured /^+ as (280 ± 19 ± 28 ± 34) MeV [22], and we are consistent with these predictions as 
well. We look forward to more data to improve the precision. 
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